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COMPARATIVE CRYSTAL FIELD ANALYSIS OF THE Ni** 
ENERGY LEVEL SCHEME IN CdCh, CdBr2, Cdlz CRYSTALS 


Mikhail G. BRIK!, Nicolae M. AVRAM?, Calin N. AVRAM? 


Abstract. The exchange charge model of crystal field was used to calculate the crystal field 
splitting of all Ni?* energy levels (3d° electron configuration) in three structurally similar 
crystals of CdCl2, CdBr2, and Cdl. The influence of the nature of the ligands and interionic 
distances on the values of the crystal field parameters and location of the split energy levels 
and relative intervals between them was identified. Comparison of the calculated energy 
levels with the experimental results available in the literature yielded good agreement. 
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1. Introduction 


The transition metal ions with the unfilled 3d electron shell are widely used in 
various applications, such as lasing [1,2], lighting [3], optical thermometry [4,5] 
etc. Their absorption and emission spectra are essentially dependent on the host 
matrix, since the electronic states arising from the open 3d electron shell strongly 
interact with the nearest crystalline environment. 


The theory of crystal field is a powerful tool that has been successfully applied to 
analyze the optical spectra of the transition metal and rare earth ions in crystals [6]. 
Different models of crystal field that have been developed so far allow to establish 
firm and deep connections between the crystal structure (the nearest surrounding of 
an impurity ion) and chemical nature of the surrounding ions on the one hand and 
the overall splitting of the impurity ion’s energy levels on the other hand. 


In this connection, comparative studies of the crystal field effects and impurity ions 
energy levels splitting for a number of similar or different systems are of special 
importance. For example, calculations of energy level schemes of the same impurity 
ion in several isostructural crystals highlight the role of the host material in the 
formation of the impurity ion’s spectra [7, 8]. Alternatively, calculations of spectra 
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of several different impurity ions in the same host help to identify the whole spectral 
range of absorption and emission that can be realized for a given host [9, 10]. 


In the present paper we report on the results of applications of the exchange charge 
model of crystal field [11] to the calculations of the Ni’* energy levels in three 
structurally close crystals — CdClo, CdBr2, and CdI2. Comparative discussion of the 
crystal field parameters and splitting of the degenerated Ni’* energy levels is given. 
A clear trend in variation of the above-mentioned parameters and the role played 
by the type of ligands in these detected changes is identified and explained. 


2. Crystal structure of CdChk, CdBr2, CdI2 and method of calculations 


CdCl2 and CdBrz crystallize in the R-3m space group (No. 166) and Cdl — in the 
P63mc space group (No. 186). The summary of the structural properties of these 
crystals is given in Table 1. 


Table 1) Structural parameters of the CdX2 (X = Cl, Br, I) crystals [12]. 


CdCl CdaBr2 Cdlz 
Space group No. 166 166 186 
a=b,A 3.85 3.95 4.24 
c, A 17.46 18.67 13.67 
ea 120 120 120 
Cd —X distance, A 2.6567 2.7607 2.9854 
180 (x3); 180 (x3); 180 («3); 
X—Cd-—X angles, ° 92.87 (x6); 91.35 (x6); 90.49 (x6); 
87.13 (x6) 88.65 (x6) 89.51 (<6) 


It is seen that the Cd — ligand distance increases in the Cl — Br—I direction. In each 
of these crystals the Cd ions are six-fold coordinated at the centers of a slightly 
deformed octahedron. The deformation of the coordination octahedra can be 
assessed in terms of the interionic distances and internal angles in the octahedron. 
The six Cd — X distances in each of these halides are equal, but there exists an 
angular distortion. In the ideal octahedron there are three X — Cd — X angles equal 
to 180 ° and twelve X — Cd — X angles equal to 90 °. In the considered halides, 
however, the angles which should be 90 ° deviate from this value; this deviation is 
the greatest in CdCl» and the smallest in Cdly. The local symmetry at the Cd”* site 
is Dsa. 


Fig. 1 shows the unit cells of CdBr2 (which is isostructural to that one of CdClz) 
and CdIz. A common feature of all three halides is that they have a well-pronounced 
layered structure; the layers made of the CdX6 octahedra are perpendicular to the c 
crystallographic axis. After doping, the Ni’* ions occupy the Cd** sites without 
charge compensation. 
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Fig. 1. a) One unit cell of CdBr2 (a) and Cdl2 (b). Coordination octahedra around Zr ions are shown. 
Drawn with VESTA [13]. Octahedral coordination around Cd ions is shown. 


The Ni** energy levels were calculated as the eigenvalues of the crystal field 
Hamiltonian 


H= Dneaia a By Op (1) 


Here OF are the linear combinations of the irreducible tensor operators acting on 
the angular parts of the impurity ion’s wave functions (exact definition of the 
operators used in the exchange charge model (ECM) can be found in Refs. [11, 14), 
and Bkare the crystal field parameters (CFPs), which can be calculated from the 
experimentally obtained crystallographic data. The Hamiltonian (1) should be 
diagonalized in the space spanned by all wave functions of the free ion’s LS terms 
arising from the Coulomb interaction between electrons of an impurity ion. In the 
ECM framework the CFPs are defined as a sum of two contributions [11]: 


k — pk k 
By = Bua + By. 
where 


vi Oud 
Big = —Kye(r?) igi ape (3) 


2(2pt1 V5 (0 ,0i) 
Bys = Ke? SRO Y(G.5(s);7 + GyS(0)i? + YpGnS(m)°) (a) 


The first term Be 6 describes the electrostatic interaction between the central 
impurity ion and the crystal lattice ions enumerated by an index 7 with electrical 
charges qi (in the units of the proton charge) and spherical coordinates R;, 6;, p; (the 
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reference system is centered at the impurity ion itself). The averaged values(r”) (r 
is the radial coordinate of the d electrons from the unfilled electron shell of an 
impurity ion) can be either found in the literature or calculated numerically, using 
the radial parts of the corresponding ion’s wave functions. The numerical 
factors Ke Yp, the explicit expressions for the polynomials ve and the definitions 
of the operators of can all be found in Ref. [11]. 


The second term of Eq. (2) Bes is proportional to the overlap between the wave 
functions of the central ion and ligands. The S(s),S5(o),S(a) terms correspond to 
the overlap integrals between the d-functions of the central ion and p-, s-functions 
of the ligands: S(s) = (d0|s0),S(a) = (d0|p0),S() = (d1|p1). The G,, Go, G; 
entries are the ECM dimensionless adjustable parameters. Their numerical values 
are determined from the positions of the first three absorption bands in the 
experimental spectrum. A good approximation is to assume them to be equal, G,; = 
G, = G, = G; then only the lowest in energy absorption band is sufficient to find 
the value of G [11]. The summation in Eq. (4) includes the nearest neighbors of an 
impurity ion only (i.e. six/four ligands in the case of an octahedral/tetrahedral 
impurity center), because the overlap of the wave functions decreases rapidly with 
the distance and is important for the nearest neighbors only. 


Despite a small number of fitting parameters employed in the ECM, it is possible 
to calculate the crystal field parameters and energy levels of impurities in crystals 
without invoking any assumptions about the impurity center symmetry, and with 
high accuracy. Many examples of successful applications of the ECM to the 
calculations of energy level of the rare earth ions [11, '* and references therein] and 
transition metal ions [Eroare! Marcaj in document nedefinit.,15,16, 17 and 
references therein] can be found in the literature. 


3. Ni** ions: energy levels in a free state and in a crystal fields of On and D3a 
symmetries 


Ni** ions have the following electron configuration: 1s?2s?2p°3s73p°3d*. The 
completely filled electron shell are not optically active, and the main optical 
properties are determined by the 3d° electron configuration (if the high-lying 4s, 4p 
states are not considered). 


The Coulomb repulsion between the 3d electrons gives rise to 5 LS terms: two spin- 
triplets 3F (the ground state) and °P and three spin-singlets 'D, 'G, !S. Their relative 
energies are expressed in terms of the Racah parameters for a free ion Bo, Co as 
follows (if the energy of the °F term is taken as zero): 15Bo for 37P, 5Bo + 2Co for 
'D, 12Bo + 2Co for 'G, and 22Bo + 7Co for |S (the subscript “0” indicates that these 
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values are for free ions; when the ions with unfilled electron shells are embedded 
into the solids, the values of the Racah parameter are reduced due to the 
nephelauxetic effect [18]). 


In a crystal field these energy levels are split into a number of sublevels, whose 
degree of degeneracy and symmetry properties of their wave functions can be 
analyzed with the help of the group theory [1, 17]. Table 2 below shows the splitting 
of these LS terms for the On group and Dsa group (the latter one is the point 
symmetry of the Ni** ions in the considered hosts). 


Table 2) Splitting of the orbitally degenerated states (the values of the orbital momenta L are given 
in the parenthesis) of the 3d° electron configuration in the crystal fields with the Opn and D3q local 
symmetries. The entries in the parenthesis in the D3qa column show the splitting of the corresponding 
entries in the O, column. 


On D3a 
S (L=0) Aig Al 
P (L=1) Tie A2t+E 
D (L=1) Toe + Eg (Ai+ E) + (E) 
F (L=3) Aog + Tig + Trg (Az) + (A2 + E) + (Ai +E) 
G (L=4) Aig t+ Ey t+ Tigt Tre (Ai) + (E) + (42 + E) + (Ai +) 


The quantitative analysis of the energy levels of d ions in crystal fields with the 
cubic (On or Ta) symmetry is based on the Tanabe-Sugano diagrams ['*]. Since three 
parameters are needed in this case to calculate all energy levels positions (the crystal 
field strength Dg, Racah parameters B and C), these diagrams express the energy 
of the crystal field states in the units of the Racah parameter B versus the Dq/B ratio 
for the fixed C/B ratio. 


Fig. 2 shows the Tanabe-Sugano diagram for the d® electron configuration in the 
octahedral crystal field for the C/B ratio equal to 4.5. When Dq/B = 0, the free ion 
energy level scheme is obtained. When Dq/B > 0, all energy levels become split in 
accordance with the data from Table 2. The energy interval between the ground 
state *A2, and the first excited spin-triplet state >72, (both come from the ground 
term °F) is equal in this case to 10Dq value [19], which gives an easy way of 
experimental determination of the crystal field strength for the d? and d® electron 
configurations. 
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E/B 


Dq/B 


Fig. 2. The Tanabe-Sugano diagram for the d® electron configuration in the octahedral crystal field. 
The solid and the dashed lines show the spin-triplet and spin-singlet states, respectively. 


4. Results of calculations 


The structural data from Ref. [12] were used to create large clusters consisting of 
25221 ions in CdCl2 and CdBr2 and 24338 ions in Cdl2, to ensure convergence of 
the crystal lattice sums needed to calculate the point charge contributions to the 
CFP. The overlap integrals between the wave functions of the Ni?* and Cl, Br, I 
ions were calculated numerically using the radial parts of their wave functions from 
Ref. [20] and were approximated by the exponential functions, whose equations are 
given in Table 3 [7]. 


The non-zero calculated CFP are collected in Table 4. The trigonal symmetry of the 
crystallographic site occupied by the Ni** ions is reflected by the structure of the 
crystal field Hamiltonian, in which only three CFP - B?,B?, B? differ from zero. 
The eigenvalues of the crystal field Hamiltonian (1) with the CFP from Table 4 are 
collected in Table 5, in comparison with the available experimental and theoretical 
data. 
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Table 3) The Ni** - X (X = Cl, Br, I) overlap integrals (R is the Ni** - X" distance in atomic units). . 


Ni?* - Cl 


Ni2* - Br 


Ni?* - I 


S(s) = (d0|s0) 


-1.82070 exp(-0.85599R) 


1.53760 exp(-0.79983R) 


1.22140 exp(-0.72022R) 


S(a) = (d0|p0) 


0.81606 exp(-0.59570 R) 


0.57804 exp(-0.51055R) 


0.49431 exp(-0.46882R) 


S(t) = (d1|p1) 


1.54200 exp(-0.85907 R) 


1.42630 exp(-0.81434R) 


1.33490 exp(-0.76824 R) 


Table 4) Calculated CFP (in cm’!, Stevens’s normalization) for Ni* ions in CdCl2, CdBr2 and Cdl 
crystals. Q, S, and Total stand for the Baas Bys and total CFP, correspondingly. The values of the 
non-dimensional ECM parameter G and Racah parameters B, C (in cm!) are also given. 


CdClz CdBr2 Cdl, 

Q S Total Q S Total Q S Total 
BS 1261 -540 721 1341 -239 1102 1152 -83 1069 
Be -72 -1096 | -1168 -63 -1071 | -1134 -45 -1070 | -1115 
B3 -1669 | -34555 | -36224 | -1360 | -31790 | -33150 | -911 | -30776 | -31687 
G 15.12 15.45 21.0 
B 800 600 825 
C 3470 3670 3250 


It is seen from Table 4 that the Racah parameters are considerably reduced to the 
free Ni** ion values (Bo = 1068 cm!, Co = 4457 cm’! [21]), which indicates high 
degree of covalence in all considered materials. 


The overall agreement between the calculated and experimentally detected energy 
levels shown in Table 5 is good. It should be emphasized that only one fitting 
parameter G of the ECM is sufficient to reach such a good correlation with the 
experimental data; no assumptions for the local symmetry of the impurity center 
was made, but only the crystal structure data were used as an input for the crystal 
field calculations. The obtained CFP values also allowed to analyze the crystal field 
splitting in terms of the D3za point group irreducible representations by using the 
symmetry descent approach from the parent On to the actual D3aq symmetry. 
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Table 5) Calculated energy levels (in cm™, Stevens’s normalization) for Ni* ions in CdCl2, CdBro 
and Cdl crystals. The orbital doublet states are denoted with an asterisk. 


CdCl CdBr2 Cdlz 
Exp.° Calc.” | Calc.¢ Exp.“ Calc.” | Calc.° Exp.4 Calc.” Calc.4 

3 

on 4 | 0 0 0 0 0 0 0 0 0 
3T¢ 3A 7200 6983", | 7140 6645 6456, | 6645 | 6350- | 6240", | 6320- 
CF) 3E 7457 7231 7014 6659 | 7100 6728 7140 
>T ig 3E 11365 11849", | 11869 | 10616, | 10770, | 10982 | 11200 | 10718, | 10520- 
CF) 3A. 12116 | 12029 | 11416 | 10887* | 11005 10884" | 11300 
1 * % 

E, 1 , 11820 | 11826 11218 | 12540 | 12514 12510 
('D) E | 12920 | 12845 12264 

Te TAY 19441", | 19005 | 16800 | 17981", | 17550 18292*, 
('D) 'E 19923 | 19179 18271 | 18547 18899 

‘Ale ‘4, | 19200 | 20602 | 19748 18116 | 18760 | 18072 20149 
('G) 

3T 1 342 | 21700, | 21185", | 20653 | 20580 | 17608", | 19228 20276", 

CP) 3E | 22200 | 22256 | 21090 19273 | 19293 21842 

Tie "A> 23560", | 22843 21038", | 21003 22621", 
('G) IE 23997 | 22931 21554 | 21017 23128 

1 * % 
ie IE 28702" | 28232 25826" | 25990 27016 

Te IE 28713, | 28604 26017, | 26316 27019, 

('G) Ay 29885" | 28798 27419" | 26275 28486" 

1 

ne 1A, 52147 | 50313 48419 | 46108 50015 


@ Ref. [22]; ° This work; ° Ref. [23] 4 Ref. [24] 


The calculated values of Dg (which were estimated from the barycenters of the 
states coming from the *7>» state) are 714 cm’! for CdCl:Ni**, 664 cm’! for 
CdBr:Ni**, and 640 cm’! for Cdlz:Ni?*. These values are in good agreement with 
the experimental estimations of 720 cm’! for CdCh:Ni** [22], 640 cm?! for 
CdBr2:Ni** [22] and 666 cm! for Cdly:Ni** [24]. 


The decrease of the calculated Dg parameter in the CdCl2 — CdBr2 — CdI2 sequence 
follows the increase of the interionic distances in the same direction. 


The Dq/B ratios are equal to 0.89 for CdClo:Ni”*, 1.10 for CdBr2:Ni**, and 0.78 for 
Cdln:Ni**, which indicates that a weak crystal field case is realized in each of the 
considered systems. 


COMPARATIVE CRYSTAL FIELD ANALYSIS OF THE Ni?* ENERGY LEVEL SCHEME IN 
CdCl, CdBr2, CdIlz CRYSTALS 27 


5. Conclusions 


Three cadmium halide crystals — CdCl2, CdBr2 and Cdl2z, doped with the Ni** ions 
— were considered in this paper. The crystal field parameters for the Ni?* ions were 
calculated using the exchange charge model of crystal field, and the obtained values 
were used to diagonalize the crystal field Hamiltonian to get the Ni** energy levels. 
The D3a local symmetry of the crystallographic site occupied by the impurity ions 
in all three studied compounds is reflected in the structure of the crystal field 
Hamiltonian and the pattern of the energy levels splitting. The agreement between 
the calculated and experimentally detected energy levels for these systems was 
good. A decrease of the crystal field strength with increased anion atomic number 
was explained by increased interionic distances. 
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